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deduce the Statistical Distance metric. By computing the curvature of this 
metric we can identify regions of parameter space most sensitive to changes 
in these parameters and thus lead to optimum detection statistics. 

I. INTRODUCTION 



There has been an expanding interest recently in the geometry of quantum state space 
|jl|||. Although some general features of the geometry are known, little has been learned 
concerning the details of the geometry of high dimensional pure and impure states. The ini- 
tial discovery of a geometric phase by Berry M was interpreted by Simon || as the holonomy 



transformation in parallel transporting the adiabatic eigenstate in parameter space || . Since 
then the concept of the geometric phase has been broadened to cope with non-adiabatic, 
non-cyclic and non- unitary evolutions 0. However, although the formal understanding of 
the geometric phase has progressed, knowledge of the underlying geometry described by this 
phase has not. This is mainly due to computational difficulties in calculating metric tensors. 
These metric tensors are functionals of infinite dimensional density matrix operators. In 
this paper we calculate the metric and curvature of the parameter space of squeezed ther- 
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mal quantum states. Using a basic understanding of the quantum metric from a statistical 
inference viewpoint we identify regions of parameter space which yield large changes in the 
quantum state making its determination easier in such parameter regime. In section II we 
review the work done to date in uncovering the geometry of a quantum state and introduce 
the Bures metric. The main results of this paper is the method of calculation and final 
results in section III. We have tried to keep this short as short as possible while including 
enough detail for the reader to reproduce their own calculation. 

II. REVIEW 

In this section we review the basics of the geometry of quantum states and introduce 
concepts associated with the natural extension of the Fubini-Study metric to impure density 
matrices and in particular, the Bures metric. 

For pure quantum states the geometry is CP n and is essentially the geometry of the 
horizontal section of the fibre bundle over the space of pure states with a fibre group U(l). 
The connection defining this section is "natural" in that the resulting metric and distance 
functions are invariant under a global change in phase of the states involved. More precisely, 
the expectation value of any operator in the state is unchanged under the action \tp) — > 
e i9 \^). We can thus split up the space of pure states into conjugacy classes under the U(l) 
action and denote the class formed from the state |?/>), the ray at We can define a 
distance between two rays as the smallest transition probability between any two elements 
in the separate rays ie. D FS = inf || \ipi) — e ie \ip2)\\ 2 - Extreamising over the relative phase 
9 we obtain the well known Fubini-Study distance for pure states D FS = 2(1 — | (■?/>! |?/>2)| 2 ). 
One can show that the geometry is Kahler and the metric, ds 2 FS , is the Hessian of a suitable 
Kahler potential ||. This Riemannian metric arises in calculations of Berry's phase and 
generalisations thereof J7| and in aspects of quantum distinguishability fll2| . The geometry 
of impure quantum states has received little attention. A Riemannian metric for classical 
probability distributions was obtained independently by Wootters || and Campbell ||. A 
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transition amplitude between two impure quantum states was discovered by Bures |10|. This 



amplitude and the related metric has been studied at length with regard to geometric phase 



by Uhlmann The same metric has also been obtained in other work relating to the 



optimal statistical distinguishability between two quantum states | fL2|j . 

Although some formal work has been done on the geometry of impure quantum states 
few concrete results concerning the details of the metric have been found. This is due to 
the technical difficulties in computing the Bures (or Statistical Distance) metric. Before 
giving the formula we will outline briefly the origins of this metric following Uhlmann. The 
derivation follows the above argument for the Fubini-Study metric on pure states. Beginning 
with an impure state p one purifies this state by enlarging the Hilbert space into a Hilbert- 
Schmidt space through the "square-root" of p, i.e. H — > H ext = H®H* where p = WW*, W 
is pure in T~C ext and where Tr WW* < oo. The "square-root" W, of p is defined up to right 
multiplication by an arbitrary unitary operator V. We again have a fibre bundle structure 
with base space y/p and fibre Ry, where Ry is right multiplication by V. The natural 
distance in H ext is the Hilbert-Schmidt metric d 2 (W 1 ,W 2 ) = Tr (W x - W 2 ){W l - W 2 )*. 
This gives a natural connection on the bundle and one can again define a distance between 
two fibers to be the smallest Hilbert-Schmidt distance between elements of the fibres. The 
solution to the extreamisation of VK(A) where A is an afline parameter, is W = GW where 
G = G* and • = d/dX. The induced metric on the horizontal section is just 

(W, W)hs = {GW, GW)hs = Tr G 2 p = ^Tr Gp . (1) 

The extreamised W(X) is parallel transported with respect to the natural connection and 
gives rise to an evolution for p(A) in Tt which obeys 

p = Gp + P G . (2) 

The Bures distance which results from the extreamisation can be written as 

DU Pl ,p 2 ) = 2 



1-TrV/^W 2 



(3) 
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while the infinitesimal Riemannian metric resulting from this distance is 

ds 2 B = TrG 2 p = ^Tr Gdp , (4) 

where 

dp = Gp + pG . (5) 
This metric is also known as the Statistical Distance metric and is symmetric in pi, p2 [fL2 



Investigation into the detailed structure of the Bures distance has been hampered by the 
complicated square-root factors in (|^). The distance and metric have been calculated for 



the spin-1/2 system [|T^] and the spin-1 system p4jfi5|j . It was found that the geometry of 
spin-1/2 state space was of constant curvature. However, the geometry of the spin-1 state 
space possessed a non-constant curvature. It was further proved in that the geometry 



of state space for spin-n is not of constant curvature and not even locally symmetric for 
In + 1 > 3. To directly solve for the Riemannian metric (|J) one must solve the matrix 
Lyapunov equation (13). For quantum systems possessing a finite dimensional representation 



the method of annihilating polynomials can be used to solve the Lyapunov equation JT? . 
This becomes prohibitive for n > 3 and results in non-unique expressions with respect to the 
parametrisations chosen for S p. Other solution methods are available but are again difficult 
to compute i.e. recursive solutions. 

In the following we will first calculate the Bures distance between two undisplaced ther- 
mal squeezed states and from this derive the associated Riemannian metric. Calculating the 
curvature of this metric we find the space is not of constant curvature and can interpret this 
curvature as a measure of optimal quantum distinguishability between the states. 

III. BURES DISTANCE 



From the work of Bures and Uhlmann p4|TT| the transition amplitude between two 
quantum states may be written as 



4 



D 2 B ( Pl ,p 2 ) = 2 



1 - TrVp} /2 p 2 pl /2 



(6) 



Due to the complexities of computing the trace, studies of this transition amplitude have 
concentrated only on finite dimensional examples with concrete results for dimension 2 |14 



and 3 ||15|| . In the following we will compute the transition amplitude between two thermal 
squeezed states with density matrices parametrised in the form: 

p(f3,r,9)=ZS(r,e)T((3)tf(r,9) , (7) 

where 

S(r,9) = exp (CK + - (*K„) , 
T((3) = exp(-(3K ) , 
C = re ie , 



and 
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K+ = -a ] , K- = \a 2 , K = -(a^a + - 



[K ,K±] = ±K±, [K.,K + ] = 2K 



Here S(r, 9) is the one photon squeeze operator, a is the single mode annihilation operator, 
Z is chosen so that Tr(p) = 1, and (K ,K±) are the generators of the SU(1,1) group. 
Equation (^) thus represents an undisplaced squeezed thermal state. We have written the 
density matrix in the Schur form p = UTU^ where U is unitary and T is diagonal in the 
eigenbasis of Kq. This decomposition is relatively straightforward in the case of gaussian p 
fl8|| . However, we find that the following arguments do not seem to hold if we expand the 
states considered to include displaced thermal squeezed states. We are thus restricted to 
density matrices continuously paramterised by three variables (/3,r,9). 

In the following this Shur factorisation will play a central role. With this factorisation 
we can easily define the square root of a positive operator. A possible alternative method 



is to represent the squeeze and thermal operators as 2x2 matrix represntations of SU(1,1). 
However, it is unclear to the author at this time how one can consistently define the square 
root in this representation without again forming the Shur decomposition of the matrix 
representation. We now outline how the Shur decomposition yields the square root of a 
positive operator p. 

Through an insertion of unity, the Schur factorisation of p yields a diagonal representation 
of the state. Choosing orthogonal eigenstates \\), such that T(/3)\\i) = Aj(/3)|Aj) and 
ZJ2i Aj = 1, we can insert the resolution of unity J2i = 1 into (0) to get 

p= ZUTU ] = ZUT\U ] 
= Zj2UT\\ i )(\ l \W = ZY,kU\\i)(\i\tf 

i i 

= Z (ui\ , 

i 

where \ui) = U\Xi), (ui\uj) = Sij and Tr(p) = 1. Thus we have diagonalised p over a 
complete orthonormal set of states with corresponding probabilities Pi = Z\. It is now an 
easy matter to find p 1 / 2 : 

p 1 ' 2 = W = Y J Pl /2 \u l ){u l \V , (8) 

i 

where V is an arbitrary unitary operator and p = WW^ . Essentially, V encodes the ambi- 
guity in taking the square root of an infinite dimensional operator. To evaluate the trace in 
(§) we need only compute p X \ 2 . From the alternate definition of the Bures Distance 

D 2 B = (VpT, Vp~2)hs = inf Tr(W 1 - W 2 ){W 1 - W 2 Y , (9) 

where p\ = WiWi, p 2 = W 2 W 2 , we see that the definition (||) is invariant under the 
transformation Wi — > WiV where VV^ = 1. By right multiplying W\ and W 2 by Vi in (||) 
we can shift away the V\ dependence of yfp to get 

pi 2 = E p iK)< u <l ' ( 10 ) 

i 

Since we need not calculate p X 2 2 we can ignore the unitary V 2 Vl appearing in p 1 ^ 2 . 
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Let us now summarise the manipulations needed to calculate the complicated trace factor 
in (^). Given p\ we can now compute ^fp taking the positive section for the square root. 

1 /2 1 /2 

Using p 2 we form the operator A = p/ p 2 p/ and, using the Baker-Campbell-Hausdorff 
identities, we rearrange A into Schur form 

A = U A T A U\ . (11) 

In this representation we can easily compute the square root, again taking the positive 
section. All that remains is to take the trace. The rearrangement of A into Schur form 
is not trivial. The operators Ua and T A can be found only in the case of undisplaced 
squeezed thermal states. For displaced states a Schur resolution was not found through 
BCH disentangling. It may be that in the case of displaced states the positive section for 
the square-root is not the correct ansatz and the more general form (H) is needed. We will 
not address this here and will only consider states of the form (0) which do result in Schur 
decompositions for A. 



I 1 /2 1 /2 

We now calculate Tryp/ P2P1 where 

Pl = ZiSin, 9i) exp(-Air )5 t (r l , 9$ , (12) 

with normalisation Tr(p) = 1 or, Z = 2sinh/?j/4. Writing Tj = exp(— (3iK ) we must 
rearrange P\ 2 P2P\ 2 to have the Schur form: 

p\ ,2 p 2P \ 12 = Z 1 Z 2 S 1 Tl /2 S\S 2 T 2 StS 1 T^St 

= Z 1 Z 2 S 1 S 3 T 3 Stst = A . (13) 

Taking the positive square root and trace gives: 



TyVI = dZ x Z 2 Tr e -^ /2Ko 



sinh (3i /4 sinh (3 2 /4 



sinh (3 3 /8 

We must now use Baker-Campbell-Hausdorff identities to express (3^ in terms of 
(Pi, p 2 , r 2 , 6*i, 9 2 ). This first step is to collapse the product S\S 2 in flT3"D into a single 
squeeze operator. This is accomplished through the identity pi] 



S^n, Oi)S(r 2 , 6 2 ) = e~^S(f, 6 - <j>)R{<j>) 



(15) 



where R(4>) = e lr ^ K °, is the rotation operator while the parameters (r^, 0») are related to 
(f, 0, 0) through 



(16) 



where 



a 



r6 



cosh r e sinh r 



-2i0 gjj-Qi r cosh r 



(17) 



and (T3 is the third Pauli matrix. Collapsing this product in ([[3D gives 



(18) 



However, RT 2 R^ = T 2 since i?, T 2 commute. The factor exp(i</>) cancels since it is a scalar. 
We thus have 



P 



TpipT = z x z 2 s x t\ I2 st 2 s^t{ 12 s\ . 



(19) 



We now must rearrange the product T^ 2 ST 2 S^T^ 2 S\ into Schur form. To do this we use 
Baker-Campbell-Hausdorff disentangling. The particular method we use was outlined in 
TT| . Using the faithful group representation of SU(1,1) where 



1 









1 




K_ = 




K „ = i 

2 









-1 


-1 



K + 

we can express the operator S as 

cosh 7 e 10 sinh 7 

5 = exp(C^ + - CK-) = 

e l6, sinli7 cosh 7 
where 7 = y((* and (9 = y(/(*. The operator Tj is represented as 



(20) 



(21) 



Ti 



e -A/2 
e^ 2 



(22) 
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To re-express the product TiST 2 S^Ti as ST$S^ we represent each operator as a 2x2 matrix, 
multiply and compare the resulting entries to obtain 

cosh (3 3 /2 = 

cosh 2 7 cosh((/3i + /3 2 )/2) - sinh 2 7 cosh((/?i - (3 2 )/2) . 



Denoting 



Y = cosh 2 /? 3 /4 = 

cosh 2 7 cosh 2 + /5 2 )/4)) - sinh 2 7 cosh 2 ((/?!- /? 2 )/4) 



and inserting this back into (14) we get 



TrV/^W} 72 



2sinh/V4 sinh/3 2 /4 



Y — 1 



All that remains is to express 7 in terms of (r 1; r 2 ,9i,9 2 ). From ( O ) we obtain 



cosh 2 7 = cos 2 A9 cosh 2 Ar + sin 2 A9 cosh 2 Er 



(23) 



(24) 



where A9 = 9± — 9 2 , Ar = 77 — r 2 and Sr = 77 + r 2 . Defining /3 + = + (3 2 )/4 and 
/5_ = (A - (3 2 )/A we finally get 



and 



Y = cosh 2 £3/4 



cos 2 A9 



cosh 2 Ar cosh 2 j3 + — sinh 2 Ar cosh 2 /3_ 



+ sin 2 A# cosh 2 Sr cosh 2 /3 + — sinh 2 Sr cosh 2 



T/2 V2 _ 2sinh/?i/4sinh/3 2 /4 



TrV Pi P2P1 



Y - 1 



For the case A# = Ar = 0, ie. only a change in temperature, (^) gives 



D 2 B {p{(3 1 ),p{(3 2 )) = 2 



sinh /?i/4 sinh/9 2 /4 



sinh( 



P1+P2 



(25) 



(26) 



Equation (|2q) gives the Bures distance between two thermal states with temperatures 
proportional to and 1 / /?2- The distance function (|26| ) (or more generally (|6])) is a proper 
distance function on the space of states. However, the resulting form (|26| ) is clearly not a 
distance function arising from a local metric structure defined on the parameter space. The 
restriction of the Bures distance to pure states, the Fubini-Study distance, is derivable from 
a local metric ie. D 2 FS = 2(1 — | (^ll^) | 2 ) = 2cosA# where A6 is the angle between the 
two Hilbert space vectors I V^)- From Uhlmann's derivation of the Bures distance as 

the minimum Fubini-Study distance between the purifications of y/pi and ^J~p~2 in the larger 
Hilbert-Schmidt space we see that the Bures distance arises from a metric structure in this 
larger Hilbert-Schmidt space. To derive this metric we can proceed in two ways. We can 
use standard perturbation analysis to evaluate ds\ = D B (p, p + Sp). This was done in |TJJ 
with the result 

5D% = ds% = - £ ^Ml^ML t (27) 
2 ^ Pi + Pj 



where |itj) are the eigenstates of p. This metric also appears in |12 and is known there as 
the Statistical Distance metric. The Bures or Statistical Distance metric has been mostly 
studied in the case of pure states while for impure states only a few finite dimensional 
examples have been examined [13-T5|. Since we know the orthogonal eigenstates of p we 



could compute fl27|) explicitly. This is simply done in the case Ar = A# = but becomes 
quite tedious otherwise. Instead we note that 



t=o 



ds 2 B = g a pdx a dx 13 = 

^D B (p(P, r, 6), p(P + t5fi, r + t5r, 9 + t56)) 
Using (|25|), (HI) and (p8|) one can eventually obtain 



ds 2 B — ds 2 SD 
I 

= - [1 + sech/3/2] (dr 2 + sinh 2 (2r)^ 2 ) 



(28) 



64sinh 2 /3/4 
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d(3 2 , (29) 



where the subscripts indicate Bures or Statistical Distance. The metric may be simplified 
somewhat by defining exp(— 2u) = tanh/9/8 to give 



ds 2 B — ds 2 SD 



1 r 



1 + tanh u 



(dr 2 + sinh 2 (2r)dfl 2 ) + du 2 . (30) 



2 

As a check we have directly computed the df3 2 contribution from (|27| ) while for pure states 



(f3 — ► oo) the metric reduces to the known form pll. Equipped with the metric one 



can compute geometrical quantities such as the scalar curvature 

8 (cosh 2 /3/4 + 12 sinh 4 (3/ A) 
RsD = crf^ ■ (31) 

It is interesting to note that the scalar curvature is independent of the "unitary" parameters 

r and 9 and only depends on the "nonunitary" parameter {3. This may be understood from 

the work of Dittmann |L5|| . Dittmann shows that in finite dimensions the geometry of the 



parameter space is locally isometric to S" 1 ^ 1 x U(n)/T n . The homogeneous submanifold 
U(n)/T n , which, in our case, is parametrised by r(/3) and 6(/3), has a constant curvature 
depending only on (3. Thus, in the general case, the curvature should only be a function of 
the invariants $ ]22[ of the state p. 

IV. DISTINGUISH ABILITY MEASURE 

A physical significance can be attributed to the curvature Rsd (|31|). From a statistical 
inference viewpoint the Statistical Distance can be understood as a measure of how well one 
can, in principle, determine the parameters describing p through N arbitrary generalised 
measurements. Fro more on this viewpoint see [T21J23]]. The error 5X in estimating the 



parameter X by analysing the data obtained from N copies of p{X) is bounded by 

N{{6X)2) {^) 2 - 1 ' (32) 
where ds/dX is the rate of change of statistical Distance with respect to the parameter X 
for a single copy of p(X). Thus, if two states are separated by a Statistical Distance of ds 
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then one must perform at least N > 1/ <is 2 measurements on identically prepared copies to 
distinguish between the two. To estimate a parameter seperation between two states one 
can make use of fl32|). 

However, to go beyond distinguishing between two states and to provide a complete 
estimation of a state given a reference state one can argue that the accuracy of estimation 
of the complete state should be independant of the particular parametrisation used. The 
most natural geometric quantity which is parametrisation (coordinate) independant is the 
scalar curvature of the Statistical Distance. Following this argument, we can interpret the 
curvature Rsd as the degree of local distinguishability of the complete state i.e. if Rsd is 
large then few measurements are needed to estimate the parameters of a neighboring state 
while if Rsd is small, many measurements will be required to estimate the parameters of 
a neighboring state. We plot the behavior behavior Rsd{/3) in Figure 1. From the form of 



Rsd ([3lD we see that the degree of distinguishability diverges as (3 — > oo. Similarly, from 
(p9|) we find that the integrated Statistical Distance between a pure state and any impure 
state diverges. This feature was also seen in a similar calculation by Braunstein and Milburn 
|2"B|] . As was explained there, the degree of distinguishability, Rsd, diverges at the pure state 
boundary because only one measurement is needed to differentiate between a pure state with 
(3 = oo and an impure state with any (3 < oo. This principle was the basis of the one-shot 



clock in p3| where the parameter to be estimated was time. From we can also identify 
local maxima (at (3 = 0) and minima (at (3 = 4 cosh -1 5/\/22). The ultimate accuracy of 
parameter estimation at various values of (3 is reflected in the behavior of the degree of 
distinguishability R. This measure can serve as a guideline for optimal operating regimes in 
quantum non-demolition measurements. By targeting the measurements to operate in those 
regions of high distinguishability one can obtain information about the quantum system in 
the least number of measurements. 
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V. CONCLUSION 



In this paper we have examined the geometry of the quantum state. After reviewing the 
previous work we showed how one can extend the definition of the quantum metric to deal 
with impure quantum states. The resulting metric is known as the Bures metric. To com- 
pute the geometry of a quantum state using this Bures metric necessitates the computation 
of the square-root of a density matrix. Up till now this has only been done for very simple 
quantum systems. Using a Shur decomposition of p we have calculated the Bures distance 
and associated Riemaniann metric between two squeezed thermal states. The Riemaniann 
manifold is not of constant curvature as was suggested by spin-1/2 calculations. It was 
argued that the scalar curvature of the Riemaniann manifold can serve as a measure of the 
ultimate accuracy in determining the parameters defining a quantum state. The method 
used here can be generalised to more complicated quantum states if a Schur factorisation 
can be found. More generally, computable methods for solving the finite or infinite dimen- 
sional Lyapunov equations need to be investigated before one can understand the geometric 
structures of high dimensonal impure quantum states. 
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